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ABSTRACT 



Context. The identity of the progenitor systems of Type la supernovae (SNe la) is still uncertain. In the single-degenerate scenario, the 
interaction between the supernova blast wave and the outer layers of a main sequence companion star strips off hydrogen-rich material 
which is then mixed into the ejecta. Strong contamination of the supernova ejecta with stripped material could lead to a conflict with 
observations of SNe la. This constrains the single-degenerate progenitor model. 

Aims. In this work, our previous simulations based on simplified progenitor donor stars have been updated by adopting more realistic 
progenitor-system models that result from fully detailed, state-of-the-art binary evolution calculations. 

Methods. We use Eggleton's stellar evolution code including the optically thick accretion wind model and taking into account 
the possibility of the effects of accretion disk instabilities to obtain realistic models of companion stars for different progenitor 
systems. The impact of the supernova blast wave on these companion stars is followed in three-dimensional hydrodynamic simulations 
employing the smoothed particle hydrodynamics (SPH) code GADGET3. 

Results. For a suite of main sequence companions, we find that the mass of the material stripped from the companions range from 
0.11 M e to 0.18 M Q . The kick velocity delivered by the impact is between 51 km s -1 and 105 km s -1 . We find that the stripped 
mass and kick velocity depend on the ratio of the orbital separation to the radius of a companion, a f /R. They can be fitted in good 
approximation by a power law for a given companion model. However, we do not find a single power law relation holding for different 
companion models. This implies that the structure of the companion star is also important for the amount of stripped material. 
Conclusions. With more realistic companion star models than those employed in previous studies, our simulations show that the 
hydrogen masses stripped from companions are inconsistent with the best observational limits (< 0.0 1M Q ) derived from SN la nebular 
spectra. However, a rigorous forward modeling from the results of impact simulations with radiation transfer is required to reliably 
predict observable signatures of the stripped hydrogen and to conclusively assess the viability of the considered SN la progenitor 
scenario. 

Key words, stars: supernovae: general - hydrodynamics-binaries: close 



1. Introduction 

Type la supernovae (SNe la) play a fundamental role in astro- 
physics. They are one of the most powerful tools in cosmology 
due to their extreme luminosities and high homogeneity. Based 
on an empirical relation between the light curve s hape and the 
peak luminosity dPhillipsll 1 993t iPhillips et~ai]| 1 9991) . they can be 
used as accurate cosmic distance indicators. This le d to the dis- 
covery of the accelerating expansion of the Unive rse dRiess et alJ 
119981; iPerlmutter etlill 1 1 999l ILeibundgutll2008l) . However, de- 
spite recent progress on both, the theoretical and observational 
side, the nature of SN la progenitors a nd the physics of the ex- 
plosio n mechanisms are still unclear dHillebrandt & Niemeveii 
l2000h . 

It is widely accepted that SNe la arise from a mass- 
accreting white dwarf (C O WD) undergoing a t hermonuclear 
explosion (for a revie w see iHovle & Fowleill 19601: If inzi & Wolf 
Il967t lNomotol[T982l) . At present, the most general classifica- 
tion of progenitor scenarios are single-degenerate and double- 
degenerate models. In the single degenerate (SD) scenario, a 



white dwarf (WD) accretes hydrogen-rich material from its 
non-degenerate companion, where the companion star could be 
a main-sequence (MS) star (WD + MS channel), a slightly 
evolved subgiant star or a red giant (RG) star. When the mass of 
the WD approach es the Chandrasekhar mass limit, it explodes 

Hachisu et al l Il996b 



as a SN la (e.g. IWhelan & Ibenl 11973 



lHan & Podsiadlowskil2004l:IWang et alJ2010h . In the double de- 
generate (DD) scenario, two CO WDs spiral in and merge due 
to gravitational wave radiation, resulting in a single object with 
a mass a bove the Chandrasekha r limit, which th en may explode 
as SN la (llben & Tumkovlll984l:IWebbinll[T98l . 

However, neither the SD nor the DD models can yet be 
clearl y favored from observation or theory dMaoz & Mannuccl 
1201 lh . In the SD case, only a fairly narrow range in the accretion 
rate above 10~ 7 M yr~' is allowed in order to attain stable hy- 
drogen burning on the surface of the WD, avoiding a nova explo- 
sion. TMsjmfesjtjMicuJ^ nearby SN la 
rate annucc 3 12001 Im aoz & Mannuccil | 201 1[). Recent obser- 
vations have identified some DD binaries dNelemans et al.ll2.005b 
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iGeier et ai1l2007l) . However, only few DD systems have been 
found whose orbital per i od is short enough to merge in a Hubble 
time ( IGeier et al.ll2010t iRodriguez-Gil et al1l2010t) . In none of 
them, the combined mass exceeds the Chandrasekhar limit. In 
addition, the DD channel has been suggeste d to lead to an 
accre t ion-induced collapse ra ther than an SN la dNomoto & Iben 
11985b fS aio & Nomotolll998l) . Recently, however. iPakmor et alJ 
d2010l 1201 ll) found that violent mergers of equal-mass white 
dwarfs with masses ~ 0.9 M Q can directly trigger thermonuclear 
explosions which resemble sub-luminous 1991bg-like SNe la. 
The violent merger of two CO WDs with masses of 0.9 M Q and 
1.1 M can lead to events that reproduce obser vational charac- 
teristics of normal SNe la (IPakmor et al ] l2012bl) . 

With detailed binary population synthesis (BPS) calcula- 
tions, some research groups investigated the delay-time distri- 
butions (DTD) and birthrates of SNe la for different forma- 
tion channels. Theoretical predictions are compared with ob- 
servations to constrain th e progenitor systems of SNe la and 
their forming scenarios (|Ruiter et aO 
Wang et all [20091: IWang & Hani 120101: 



2009t iMannuccil 



Maoz & Badenes 



2009; 



2010; 



Maoz et al.l 120101) . However, there is no evidence that the SD 



scenario alone explains the shape of the observed delay-time dis- 
tribution, while this may be p ossible for the DD scen ario or a 
combination of both channels (iMennekens et al.ll2010l) . 

The obvious difference between the SD and DD scenario is 
that the companion stars in the SD channel will s urvive and may 
affect the observational di splay of the explosion dMarietta et al.l 
l2000HPakmoretal.ll2008h . In contrast, no stellar remnant exists 
in the case of the merger of two WDs in the DD channel. It is a 
promising approach to constrain the nature of SN la progenitors 
by directly searching for the surviving companion stars in rem- 
nants of SNe la. A promi nent example is Tycho Brah e's super- 
nova (SN 1572) for which lRutiz-Lapuente et al.l d2004l) presented 
a survey of the central region of its remnant, around the position 
of the explosion. They identified a subgiant star (Tycho G), sim- 
ilar to the Sun in surface temperature but with very low surface 
gravity, which moves at more than three times the mean velocity 
of the stars in the neighborhood and suggested that this might 
be the surviving compa nion of SN 1572. However, the study of 
iKerzendorf et al.ld2009l) casts some doubts on this identification. 
Recently. ISchaefer & Pagnottal (1201 2|) analyzed the center of the 
LMC SN la remnant SNR 0509-67.5 which does not contain any 
stars down to the observational limit. Thus, all possible com- 
panion star types except white dwarfs can be excluded for this 
object. 

In the SD scenario, hydrogen-rich or helium-rich circum- 
stellar material (CSM) is expected to exist around SNe la as the 
result of mass transfer from the compa nion, as well as WD winds 
dNomotolll982t lHachisu et al.lll999bh . Therefore, another indi- 
rect way of identifying SN la progenitor systems is to search for 
the material trans ferred to the accreting w hite dwarf in the C SM 
dPatat et al.l2007b . Following this approach lPatat et al.ld2007b re- 
ported detection of such CSM. They suggested that the SD sys- 
tem consist ing of a WD and RG is the favored progenitor for 
SN 2006X dPatat et al.ll2007T) . Moreover. Isternberg et al.l (1201 lh 
showed that the velocity structure of absorbing material along 
the line of sight to 35 SNe la tends to be blueshifted. Thus, they 
concluded that many SNe la in nearby spir al galaxies may orig- 
inate in SD systems dSternberg et al.ll201 lh . Note, however, that 
abundant CSM is in conflict with the missing radio signal for 
other events, e.g. SN 201 lfe dChomiuk et al.ll2012l:lHoresh et al.l 
l2"0l2h . 

SNe la are characterized by the lack of hydrogen in their 
spectra. This is explained naturally in the DD progenitor sce- 



nario. However, in the SD scenario, the companion star is typi- 
cally hydrogen-rich. Thus, the impact of a SN la explosion will 
strip off hydrogen-rich material from the companion and mix it 
into the ejecta, which may be in conflict with observations if the 
amount of hydrogen is large enough to be observable. Thus, we 
may be able to identify the progenitors of SNe la based on the 
exact amount of hydrogen stripped from the surface of the com- 
panion star. So far, there is no direct detection of the stripped hy- 
drogen. Based on high-resolution spectrosco py of the two ex tra- 
galactic SNe la SN 2005am and SN 2005cf. lLeonardld2007h es- 
ti mated an uppe r limit of 0.01 M applying the model described 
bv lMattila etal] d2005l). 

iMarietta etal] d2000t) (hereafter MOO) explored the impact 
of SN la ejecta on a variety of binary companions (MS, sub- 
giants, RGs) in the SD formation channel with two dimensional 
Eulerian hydrodynamics simulations. They found that the super- 
nova explosion can strip 0.15 M Q to 0. 17 M of material from the 
surface of MS and subgiant companions, while RGs lose almost 
their entire envelope in the impact. However, they ignored the 
effect of the mass transfer phase on the structure of the compan- 
ion star when they constructed their initial model. To investigate 
how the m ass transfer change s the mass stripping by supernova 
explosions. iMeng et alJ d2007l) used an analytical method to ap- 
proximate the mass loss in the impact. They found a lower limit 
of 0.035 M for the stripped mass, but their analytic method was 
based on an oversimplified description of th e interaction physics. 

An updated study has been presented bv lPakmor et alJ d2008l 
hereafter PRWH08) based on three-dimensional (3D) smoothed 
particle hy drodynamics (SPH) sim ulations. By mimicking the 
models of llvanova & Taariil d2004l) in their setups, they found 
stripped masses in the range from 0.01 M to 0.06 M , which is 
very close to the observational limits of ~ 0.01 M & . Therefore, 
they concluded that the SD scenario remains a valid possibil- 
ity for SN la progenitors. Although mass-loss from the compan- 
ion star was included in their study, it was modeled by removal 
of material from a main sequence star with a constant rate in 
single-star evolution code. This is an oversimplification and we 
therefore aim at reexamining the impact of supernova ejecta on 
a companion star that was modeled consistently in a detailed bi- 
nary evolution calculation. 

Very recently, Pan et al. (2012) studied the impact of SN la 
ejecta on binary companions in the SD scenario with the 
Eulerian hydrodynamics code FLASH for MS, RG and He star 
companions. They were able to quantify the amount of contami- 
nation with explosion ashes on the companion star by the super- 
nova ejecta in their simulations which might help to identify a 
companion star even a long time after the explo s ion. 

Based on the prescription of lHachisu et al.l dl999bl) for the 
mass growth of CO WDs, and including the possibility of the 
instability of an accretion disk around the WD on the evolution 
of binary systems, detailed binary evolution calcul ations have 
been p erformed for about 2400 close WD binaries bv lWang et al.l 
d2010L hereafter WLH10). They confirmed that WDs in the 
WD + MS channel with an initial mass as low as 0.61 M Q can 
accrete efficiently and reach the Chandrasekhar mass limit. Their 
calculations also showed that the disk instability could substan- 
tially increase the mass-accumulation efficiency for accreting 
WDs and cause SNe la to occur in systems with < 2 M donor 
stars. They found that the Galactic SN la birth rate from the 
WD + MS channel is ~ 1.8 x 10~ 3 yr~' according to their stan- 
dard model, which can account for 2/3 of the observed SNe la. 

In this work, we use the same method as WLH10 to carry out 
consistent binary evolution calculations for the single degener- 
ate MS channel of SNe la. With these more realistic companion 
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Table 1. Main-sequence companion star models. 



Model 1 




M Cii 


M c , f 


Pt 


a t 


a f /R 2 


'"stripped 


Vkick 


SNe la 1 




[M e ] 


[M ] 


[M B ] 


Id] 


[10"cm] 




[M e ] 


[km s- 1 ] 




ms_20a 


0.70 


2.00 


0.74 


0.98 


3.72 


3.06 


0.181 


51.01 


disk instability 


ms_20b 


0.90 


2.00 


1.17 


0.46 


2.40 


2.75 


0.105 


58.78 


disk instability 


ms_22a 


0.80 


2.20 


1.21 


0.29 


1.77 


2.72 


0.173 


105.29 


weak H-shell flash 


ms_20c 


0.80 


2.00 


1.22 


0.56 


2.74 


2.71 


0.171 


64.32 


weak H-shell flash 


ms_24a 


0.90 


2.40 


1.40 


0.33 


1.95 


2.63 


0.172 


94.95 


stable H burning 


ms_20d 


1.00 


2.00 


1.40 


0.57 


2.83 


2.63 


0.113 


53.32 


weak H-shell flash 


ms_23a 


0.90 


2.30 


1.50 


0.35 


2.08 


2.59 


0.162 


85.66 


stable H burning 


ms_24b 


1.00 


2.40 


1.88 


0.37 


2.25 


2.46 


0.116 


66.91 


stable H burning 


ms_28a 


1.10 


2.80 


2.00 


0.33 


2.11 


2.43 


0.159 


84.50 


optically thick wind 


ms_30a 


1.20 


3.00 


2.45 


0.44 


2.64 


2.33 


0.141 


65.34 


optically thick wind 



Notes. M wc i and M C j are the initial masses of WD and its companion at the beginning of mass transfer, respectively. M c _f, Pf, cif and R2 denote 
the final mass of companion star, the orbital period, the binary separation and the radius of secondary at the moment of supernova explosion, 
respectively. M stl i ppe d and Vkick are the stripped mass and the kick velocity caused by supernova impact. 

(1) All models have been named with the same way. For example, ms_20(a,b,c,d), the "ms" corresponds to CO WD + MS system, the middle 
number "20" means the M z j is 2.0 M s . The final alphabet "a", "b", "c" or "d" denote the different models with the different M c j but the same M C j. 

(2) The WD explodes as an SN la in the disk instability phase, in the optically thick wind phase, in the stable H-shell burning phase and in the weak 
H-shell flash phase, respectively. 



models, we expand and update the 3D hydrodynamical simu- 
lations performed in PRWH08 to investigate the interaction of 
SN la ejecta with MS companion stars. We then explore how 
the ejecta structure, the stripped mass and the kick velocities of 
the surviving companion depend on parameters of the progenitor 
model. Section|2]describes the codes and initial models used in 
this paper. The SPH impact simulations and numerical results for 
ten consistent MS companion models are presented in Section[3] 
The comparisons with PRWH08 and some implications of our 
simulations are discussed in Section [4] We summarize and con- 
clude in Section|5] 



2. Numerical methods and models 



2.1. Numerical codes 



We use Eggleton's stellar evolution code (lEggletonlll97lLll972t 
1 1973b to follow the detailed binary evolution of SD progeni- 
tor system^The lates t input physics have been implemented 
by IWang et all (120091 l2010h . Roche l obe overflow (RLOF) 
is trea ted in the code as described by lHan & Podsiadlowskil 
d2004l). The opa c ity tab les in our code have be e n com piled 
bv IChen & Tout! (l2007h from llglesias & Rogers! (1 19961) and 
lAlexander & Ferguson ( 19941) . We use a typical Population I 
composition with hydrogen abundance X = 0.70, helium abun- 
dance Y = 0.28 and metallicity Z = 0.02. We set the ratio 
of the typical mixing length to the local pressure scale height, 
a = I /H p, to 2, and the convective overshootin g parameter, 6 ov , 
to 0.12 dPols et al.lll997t ISchroder et al.lll997l) . which roughly 
corresponds to an overshooting length of ~ 0.25 pressure scale 
heights (H P ). 

In this paper, we start our binary evolution calculations from 
WD+MS binary systems. In a WD+MS binary system, mass 
transfer occurs through RLOF when the companion star fills its 
Roche lobe. If the mass transfer is dynamically stable, the trans- 
ferred material will form an accretion disk surrounding the WD. 
This accretion disk can become thermally unstable when the 
effective temperature in the d isk falls below the hydrogen ion - 
ization temperature ~ 6500K (Ivan Paradiisl[T996t lLasotall200lh . 
Here, as a possibility, the effect of the instability of the accretion 
disk on the evolution of WD + MS binaries has been included 
to constrain the accretion rate of the WD. If the mass-transfer 



rate, M2 , is higher than the critical mass-transfer rate for a sta- 
ble accretion disk, M cr ,disk, we will assume that the WD accretes 
smoothly at a rate M acc = \M<A. Otherwise, the disk is unstable 
and the mass-accretion rate of the WD is M acc = ikW Id, where 
d is the duty cycle, set to 0.1 in this work. 

We do not solve the stellar structure equations for the 
WD star when we construct the structure of the companion 
star for our s i mulatio n. Instead, we adopt the prescription of 
lHachisu et all (l!999bl) for the mass growth of a CO WD by ac- 
cretion of hydrogen-rich material from a companion. If the mass- 
accretion rate of the WD, M acc , is above a critical value, M C r,WD, 
we assume that hydrogen burns steadily on the surface of the 
WD and that the hydrogen-rich material is converted into he- 
lium at a rate M a -,WD, while the unprocessed matter is assumed 
to be lost from the system as an optically thick wind at a mass- 
lo ss rate M w ind = \ M2 \ ~ Mct/we )- The critical mass-accretion rate 
is (lHan & Podsiadlowskill2004 

M^wd = 5.3 x 10' 7 (L7 ^ X) (M WD /M o - 0.4) M^r" 1 . (1) 

where X is the hyd rogen mass fraction and Mwd is the mass of 
the accreting WD. iKatol (1201 ll) pointed out that this WD wind 
case may correspond to the quasi-regular transient supersoft X- 
ray source (SSS) such as V Sge. 

When |M acc | is smaller than M cl -wd, the following assump- 
tions have been adopted: 

1. If ±M cr . wd ^ \M RCC \ < M ct> wd, it is assumed that there 
is no mass loss and that hydrogen-shell burning is steady. 
In this case, before the supernova ex plosion, the system 
may be observed as the persiste nt SSS (lHachisu et al.ll2008l : 
iMeng & Yang|l2010tlKatoll20Tl . 



2. If ±M n 



*ci-,wd < \M R cc\ < jM cr y/u, hydrogen-shell burning is 
unstable, triggering very weak shell flashes, where we as- 
sume that the processed mass can be retained. Before the 
supernova explosion, thi s case may be observed as recur- 
rent nova of U Sc o-type (lHachisu et al.ll2008r.lMeng & Yang! 
l2010tlKatoll20Tl . 
3. If|M acc | < |M c1i wd, hydrogen-shell flashes are so strong that 
no mass can be accumulated by the WD. 
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Fig. 1. Time evolution of density distribution in the impact simulation with model ms_22a. For visualization (not for the simulation, 
though), we use cylindrical coordinates. The radial coordinate is averaged over angle accounting for the intrinsic symmetry of the 
star. The logarithm of density is color-coded. 



These three cases are accounted for in constructing the setup 
MS companion stars for our impact simulations. Furthermore, 
the mass-growth rate of th e WD star was linearly interpolated 
from a grid computed by iKato & Hachisul d2004h . The input 
physics in our bi nary evolution calculations is consistent with 
WLH10 (see also lHan & Podsiadlowskill2004l) . From our one- 
dimensional binary evolution calculations, we selected ten pro- 
genitor systems with a representative range of orbital periods and 
initial companion star masses. All resulting models are summa- 
rized in Table Q] 



For our hydrodynamical simulation of the impact of SN la 
ejecta on their companion stars, we use the GADGET3 code 
(ISpringel et al.ll200U ISpringelll2005l) . Originally, the GADGET 
code was intended for cosmological simulations, but it has been 
modified to make it a pplicable to stellar astrophysics problems 
dPakmor et al ] l2012ah . By using the initial parameters of the 



HCVQ scenario of MOO, PRWH08 showed that the SPH-based 
approach is capable of reproducing previous results obtained 
with a grid-based 2D scheme by MOO. This confirms that the 
SPH approach with the GADGET code captures the main dy- 
namical effects of the supernova impact on its companion star. 

2.2. Basic setup 

In our simulation, we use the same method as PRWH08 to map 
the one-dimensional profiles of density, internal energy, and nu- 
clear composition of the companion star as obtained from a bi- 
nary evolution calculation to a particle distribution suitable for 
the SPH code. Here, the smoothing length is chosen such that a 



1 In the HCV scenario, a CO WD accretes hydrogen by RLOF from 
a lower mass MS secondary. Such a system is formed when a CO WD 
is left in a close binary orbit by an earlier episode of common envelop e 
evolution in its asymptotic giant branch phase (see Marietta et al. 2000 ). 



4 



Z. W. Liu et al.: 3D simulations of interaction between SN la ejecta and their MS companions 



Table 2. Resolution test for ms_22a model. 



''star 


A/ 

"total 


"*particle 




M 

lv± stripped 








\Mn\ 






[km S" 1 ] 


50 000 


107 092 


2.41x10" 


5 


0.329 


138.17 


100 000 


214 170 


1.21x10" 


5 


0.215 


123.92 


1 000 000 


2 141 336 


1.21x10" 


6 


0.190 


111.75 


2 000 000 


4 282 671 


6.04x10 


-7 


0.185 


108.90 


4 000 000 


8 565 284 


3.02x10" 


-7 


0.175 


106.54 


6 000 000 


12 847 824 


2.01x10" 


-7 


0.173 


105.29 


10 000 000 


21 413 009 


1.21x10" 


-7 


0.173 


105.31 



Notes. A'star and iVtotai are the number of particles used to represent the 
companion star and the binary system, respectively. All particles have 
the same mass m palt i c i c . 



sphere of its radius enclosed 60 neighboring particles. The rest 
of the basic setup corresponds to that of PRWH08. 

To reduce numerical noise introduced by the mapping to en- 
sure that they are in hydrostatic equilibrium before we start the 
actual simulation, the MS companion stars are relaxed for 1 .0 x 
10 4 s (several dynamical timescales). If the relaxation succeeds, 
the velocities of the particles stay close to zero. Otherwise, we 
reject the SPH model, and redo the relaxatio n after adjusting the 
relaxation parameters dPakmor et al.ll2012ah . 

The supernova ex plosion is represented by the W7 model of 
iNomoto et al.l (Il984l) . This model has been shown to provide a 
good fit to the observational light curves of SNe la. Its total ex- 
plosion energy is 1.23 X 10 51 erg, and the average velocity of the 
ejecta is 10 4 km s . We place the supernova at a distance to the 
companion star given by the last separation of the binary sys- 
tem. The impact of the SN la ejecta on their binary companions 
is then simulated for 5000 s, at which point the mass stripped off 
from the companion star and its kick velocity due to the impact 
have reached constant values (see Section [3~!2l . 

In order to check the effect of the gravitational field of the 
WD, we run the impact simulation for the ms_22a model, in- 
cluding a 1 .4 M WD during relaxation of the companion star. 
The ms_22a model has the smallest separation (see Table[TJ) and 
should therefore be influenced most strongly by the WD. The 
results show that the companion star is basically not distorted 
due to the tidal force. The distortion should only be at the per- 
cent level in radius anyway, and our spatial resolution in the very 
outer layers of the star is not sufficient to resolve this. However, 
the mass in these outer layers is orders of magnitudes smaller 
than the total stripped mass. Therefore, we run all other simula- 
tions in this work ignoring the effect of gravitational field of the 
WD when the companion stars are relaxed. 

3. Simulations 

In this section, we present the numerical results of our SPH im- 
pact simulations for updated MS companion star models. To en- 
sure the reliability of the results from our simulations, we also 
perform a numerical convergence test for one selected compan- 
ion star model. We then explore the dependence of stripped mass 
and kick velocity on the ratios of initial binary orbital separation 
to companion radius, cif/Ro, for a given companion star. 

3. 1 . Typical evolution 

We discuss the evolution of the model ms_22a (Table[T|) as a typ- 
ical case. Figure [TJ shows the snapshots from our SPH impact 
simulation. We use a total of 6 x 10 6 SPH particles to represent 



Table 3. Fitting parameters of equation (f2| and (0. 



Models 


a 


V 


a i 


M 


ms_22a 


4.092 


-3.137 


3.84x10' 


-1.309 


ms_24a 


3.786 


-3.156 


4.11X10 7 


-1.509 


ms_28a 


2.503 


-3.095 


3.38xl0 7 


-1.573 


rp3_20a 


6.105 


-3.489 


6.05xl0 7 


-1.450 



0.40 



0.35- 
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Fig. 2. Temporal evolution of the mass loss from the compan- 
ion star for simulations with different numbers of SPH particles. 
Note that the number of particles gives the total of particles in 
the simulation (both supernova and companion star). 
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Fig. 3. Mass loss vs. different number of SPH particles used in 
ms_22a model simulations. 



the companion star, which corresponds to ~ 13 million total SPH 
particles being used in the simulation. The supernova is repre- 
sented by the W7 model and set up with an initial separation of 
1.77 x 10 11 cm as obtained from the binary evolution simulation. 

Figure[TJillustrates the density distribution of the companion 
star for model ms_22a from the onset of the supernova explosion 
until the star starts to relax again (5000 s). At the beginning of 
the impact simulation, the MS star is in equilibrium. The super- 
nova explodes on the right side of the companion star. After the 
explosion, the SN ejecta expand freely, approach the companion 
star from the right side, and a shock wave develops as they hit its 
surface (see first snapshot). The second snapshot shows the im- 
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Fig. 4. Density distribution of the companion star in the ms_22a model at 3000 s after supernova explosion with different resolutions. 
Density is color-coded logarithmically. 
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Fig. 5. Final stripped mass and kick velocity vs. the ratio of or- 
bital separation to the companion radius for different compan- 
ion models. The data are fitted by using the power law of equa- 
tion © or (O in this paper. The solid line corresponds to the 
rp3_20a model of PRWH08. Note that we use logarithmic coor- 
dinates here. All fitting parameters are given in Table[3] 



pact 100 s after the supernova explosion, when the shock starts 
to propagate through the companion star. At 400 s (third snap- 
shot), the shock wave reaches its center. The supernova ejecta 
flow around the companion and a hole forms in them. In the 
fourth snapshot (1500 s), the shock wave in the companion has 
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Fig. 6. Final stripped mass and kick velocity versus the ratio of 
separation to the radius of companion star, a^jRi, for different 
companion star models. The pentacles and filled triangles de- 
note the models simulated in this paper and models in PRWH08, 
respectively. 



passed the stellar core, and the supernova ejecta are mixed with 
hydrogen-rich material stripped from the companion star. The 
last two snapshots show the interaction at 3000 s and 5000 s af- 
ter SN explosion. This is the end of the phase of mass-stripping 
by the impact; the remnant of the companion star shrinks and 
relaxes to be almost spherical again. The mass stripped by the 
impact of the SN ejecta stays constant from this time onwards 
(see Fig. [2]). In Fig.Q] it can be seen that the companion star has 
moved to left by ~ 4.5 x 10 10 cm from 0.0 s to 5000 s due to the 
kick caused by the impact. 



3.2. Numerical convergence test 

To ensure reliability of our numerical results, we perform a con- 
vergence test. We use model ms_22a and carry out the sim- 
ulations for different resolutions ranging from 1.07 x 10 s to 
2.14 x 10 7 SPH particles (see Table|2|. For each resolution, the 
mass stripped from the companion star at different times after the 
supernova explosion is calculated (see Fig.|2|. The total unbound 
mass and the kick velocity obtained from the impact 5000 s af- 
ter explosion with different resolutions are listed in Table [2] The 
stripped mass is calculated by summing up the mass of all par- 
ticles that are not bound to the star any more, but were part of 



6 



Z. W. Liu et al.: 3D simulations of interaction between SN la ejecta and their MS companions 



Table 4. Stripped mass. 



Model 


M c , f 


a f /R 2 


AM' 


AM 1 


Difference 




[M G ] 




[M e ] 


[M Q ] 


[%] 


ms_20a 


0.74 


3.06 


0.18 


0.19 


5 


ms_20b 


1.17 


2.75 


0.11 


0.14 


27 


ms_22a 


1.21 


2.72 


0.17 


0.18 


6 


ms_20c 


1.22 


2.71 


0.17 


0.14 


18 


ms_24a 


1.40 


2.63 


0.17 


0.19 


12 


ms_20d 


1.40 


2.63 


0.11 


0.15 


36 


ms_23a 


1.50 


2.59 


0.16 


0.19 


19 


ms_24b 


1.88 


2.46 


0.12 


0.17 


42 


ms_28a 


2.00 


2.43 


0.16 


0.21 


32 


ms_30a 


2.45 


2.33 


0.14 


0.20 


43 



1 Numerical calculation of stripped mass in our simulations. 

2 Stripped mass analy tically estimated by using the method of 
IWheeler etalJ i!975h . Here, we directly use the density profiles of 
companion models from the ID stellar evolution calculations. 
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Fig. 7. Unbound mass vs. simulation time in ms_22a model. The 
solid line shows the total mass of all particles with a total energy 
(kinetic plus potential energy) larger than zero. The dash-dotted 
line also includes the internal energy in the sum. 



the star in the initial setup. In order to determine whether or not 
a particle is bound to the star, we calculated the sum of the ki- 
netic energy (positive) and potential energy (negative) for each 
particle. If the total energy is positive, the particle is not bound. 
Note that the center-of-mass motion of star is subtracted when 
calculating the kinetic energy for each particle. 

Figure[2]shows that the mass loss decreases when the number 
of the SPH particles used in the simulations increases (see also 
Fig. [3}. However, it is numerically well converged for more than 
about 8 million SPH particles in the simulation (this corresponds 
about 4 million particles in the companion star). The difference 
in stripped mass between simulations with 8.56 x 10 6 (the short 
dashed line) and 1.28 x 10 7 (the dotted line) SPH particles is 
smaller than 2%, the difference between 1 .28x 10 6 and 2. 14x 10 7 
(the narrow dash-dotted line) SPH particles is less than 1 % (see 
Table |2J. After 3000 s, the amount of unbound mass becomes 
constant in good approximation. The density distribution of the 
companion star corresponding to different resolutions is shown 
in Fig. [4] There are also no morphological differences among the 
simulations with 8.56 x 10 6 , 1.28 x 10 7 and 2.14 x 10 7 SPH par- 
ticles (the last three snapshots in Fig.|4|. Therefore, we conclude 
that it is sufficient to represent the companion stars with about 6 
million SPH particles in our SPH impact simulations. 
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Fig. 8. Initial density profiles of five companion star models at 
the onset of the SN la explosion. Model rp3_20a (dashed line) 
and rp3^0b (solid line) are from PRWH08. 



3.3. Results 

Based on a selection of ten realistic MS companion star models 
we simulate the interaction of supernova ejecta with their binary 
companions. The initial parameters of all binary systems used in 
this paper are listed in Table [T] 

The amount of hydrogen-rich material stripped from the sur- 
face of the companion stars, M stl i ppe d, and the kick velocities 
caused by the supernova impact, Vkkk, are also shown in Table [TJ 
(see also Fig.[5j. The stripped masses range from 0.1 1 to 0.18 M 
and we measure kick velocities between 5 1 and 105 kms~'. Note 
that this kick velocity is defined as the center of mass velocity of 
all particles bound to the companion star 5000 s after the explo- 
sion of the supernova. 

The largest stripped mass of our simulations is found for 
model ms_20a (0.18 M©). Further analysis shows that the hydro- 
gen in the center of this companion star has already been mostly 
depleted at the time the white dwarf explodes in a supernova, 
and its outer layers already begin to expand. This moves the star 



on its evolutionary track in the Hertzsprung-Russell (H-R) dia- 
gram towards the giant phase. The density profile of companion 
star ms_20a is shown in next section. Compared with other main 
sequence stars in our sample, this star has a very large radius and 
a higher density in the core. The shrinking of its inner core and 
the subsequent expansion of the outer layers make its envelope 
less bound. This explains why more hydrogen-rich material is 
stripped off when the supernova blast wave hits the companion 
star in this model. 

For a given companion star model, we investigate the de- 
pendence of stripped mass, M str i pP ed, and kick velocity, Vkick, on 
the ratio of orbital separation to companion star radius, af/Ro. 
All parameters but the orbital separation are kept constant (it 
means that we only change the value of the separation, «f , for the 
same companion model artificially). Figure[5]shows the stripped 
mass and kick velocity versus cif /Ri for models ms_22a, ms_24a 
and ms_28a. For a given com panion star, we find a similar re- 
lationship as previous studies ( Marietta et al.ll2000l : iMeng et ail 
120071: iPakmor et al.ll2008l: iPan et alJl2012l) . The dependence of 
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Fig. 9. Final velocity distribution of the stripped hydrogen-rich material that originally belonged to the companion star (left panel) 
and supernova ejecta (right panel) 5000 s after the supernova explosions for model ms_22a. 



the stripped mass on this ratio follows a power law to good ap- 
proximation (see Fig. [5]): 

M s tripped = flo ■ M . (2) 

Likewise, the dependence of the kick velocity, Vkkk on af \Ri can 
be fitted by a power law (see Fig. [5} 

Vkkk = a x ■ I — I cms . (3) 

Here, ao and a\ are two fitting constants, which, however, are not 
unique but depend on the companion star models. The parame- 
ters v and fi are the corresponding power law indices. Values 
for these fitting parameters are listed in Table [3] For compari- 
son, we also show the result from model rp3_20a of PRWH08. 
Comparing our models with PRWH08, all simulations follow a 
power law of some kind, but we do not obtain the same fitting 
parameters. This implies the importance of the structure of the 
companion star for the stripped mass and kick velocity in our 
impact simulations. Furthermore, we put all models used in this 
paper and the PRWH08 models together as a whole sample to 
examine the effect of Figure [6] shows the dependence of 

the stripped mass and kick velocity on the parameter, for 
all these models. However, we do not find a power law relation 
is still holding in Fig. [6] Therefore, again, it indicates the results 
of our impact simulations are also dependent on the details of 
the structure of the companion stars due to the history of mass 
transfer. 

4. Discussion 

4.1. Comparison with previous studies 

IWheeler et al.l (1 1 9751) derived a simple analytic formula for the 
calculation of stripped mass, and their results were confirmed 
by so me of the early numerical simulations (see iMarietta et al.l 
2000). Here, we also estimate the stripped mass using their ana- 
lytic method based on our MS companion models (see Table|4|. 
As shown in Table |4j the order of magnitude of stripped mass 
predicted by the analytic estimate agrees with our results, but 
is usually overestimated by the analytic formula compared to 
the result obtained from our simulations, with deviations rang- 
ing from 5% to 43%. 



Using a 1.0 M Q solar-like companion star model, 
IMarietta etal] d2000l) found about 0.15 M hydrogen-rich 
material to be stripped from the surface of MS companion by 
the impact of the supernova explosions. The mass of stripped 
hydrogen-rich material in our simulation is consistent with their 
results. Note that compared to the study of PRWH 08 our simula- 
tions b ring the results back to the original work of lMarietta et al.l 
(120001) only by chance - the setups of the companion stars 
are very different. While we construct our companion star 
models from consistent binary evolution calculations, their 
model ignores the effect of mass transfer on the structure of the 
companion star altogether. Although PRWH08 included such 
an effect in the construction of their companion stars, it was not 
done in a consistent binary evolution calculation. 

The amount of the unbound mass in our work is also consis- 
tent with the new multi-dimen sional adaptive mesh refinement 
simulations of iPan et aD d2012l) . But, again, they did not follow 
the full binary evolution but used initial conditions with a con- 
stant mass-loss rate when const ructing their MS com panion star 
models with the MESA code of iPaxton et al.l (1201 ll) . Therefore, 
the ag reement can be regarded as coincidental, too. IPan et all 
(12011 argue that the mass stripping is dominated by ablation 
in their simulations for the MS companion. They used a slightly 
different criterion for the unbound mass in which they included 
the internal energy in addition to the potential and kinetic en- 
ergy. This does not make an effect for our simulations because 
we are mainly interested in the unbound mass at late enough 
times when the internal energy is negligible to the kinetic en- 
ergy. Figure|7]shows the unbound mass of a companion star as a 
function of time after the supernova explosion in our SPH simu- 
lation. Whether or not a particle is unbound is decided by sum- 
ming its potential energy, kinetic energy and internal energy (or 
without the internal energy, which corresponds to the solid line 
in Fig. [7]i. Already 5000 s after the explosion most of the internal 
energy deposited by the impact has been converted into kinetic 
energy. Therefore, as in PRWH08, we neglect the internal en- 
ergy when we flag particles as bound or unbound at late times in 
our simulation. 

In our study we find that a minimum of 0.1 M Q of hydrogen- 
rich material being stripped from the companion star. This is sig- 
ni ficantly large r than the 0.035 M for stripped hydrogen found 
bv lMeng et atl (l2007h . This might be caused by their oversim- 
plified description of the interaction physics, e.g. the effect of 
the shock formed between the supernova ejecta and the com- 
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panion star was not calculated in their analytic model. Recently, 
iPan et al 1 (120121) argued that since they neglected the mass loss 
due to the ablation from the hot surface of the companion star, 
iMeng et all (l2007h underestimate the final unbound mass. 

PRWH08 found stripped masses in the range from 0.01 Mq 
to 0.06 M , which is very close the observational limit obtained 
bv iLeonardl (120071) . but significantly lower than the values we 
find for the models presented here. In order to determine the 
origin of this difference, we select two models, ms_20c and 
ms_22a (see Table [TJ, in comparison to the models rp3_20a and 
rp3_20b of PRWH08. The density profiles of these four mod- 
els at the moment of the explosion of the supernova are shown 
in Fig. [8] Models ms_22a and rp3_20b have the same radius of 
0.65 x 10 n cm. Although there are some small differences in 
the density profiles in the inner cores, the outer layers are nearly 
identical. In our simulations, only the properties of outer layers 
of the companions can significantly affect the results of the inter- 
action between supernova ejecta and companion star. Therefore, 
we chose these two models to carry out a comparison. PRWH08 
had set up the binary with a separation of 4.26 x 10 11 cm for 
their Model rp3_20b. They found a stripped mass of 0.01 M 
and a kick velocity of 24.1 kms~'. Our model ms_22a, how- 
ever, has an orbital separation at the time of the explosion of 
1.77 x 10 n cm only. Running this model we obtain a signif- 
icantly larger stripped mass of 0.17 M and a higher kick ve- 
locity of 105.29 km s . This is to be expected for changing the 
initial separation, as discussed above. Furthermore, we calculate 
the stripped mass and kick velocity with the same separation 
4.26 x 10 cm as for Model rp3_20b by using the power law 
relation of equation (O and (0 for our Model ms_22a. We find 
a stripped mass of 0.01 M Q and kick velocity of 33.02 km s , 
which is excellent agreement with Model rp3_20b. Therefore, 
excluding the density effect, the orbital separation at the time of 
the explosion is the primary factor to cause the difference be- 
tween the simulations of the two models ms_22a and rp3 J20b. 

The difference in orbital separations originates from different 
treatments of the progenitor evolution. The binary systems we 
examine here fill their Roche-lobe at the time of the explosion. 
This fixes the orbital separations of the binary systems. Instead 
of detailed binary evolution c alculations, PRWH08 dir ectly took 
parameters from the study of llvanova & Taam (12004 to mimic 
the effect of the binary evolution phase. Ivanova & Taaml d2004l) 
analyzed the evolution of possible SN la progenitor systems con- 
sisting of a WD and an evolved MS star. In PRWH08, all values 
adopted by mo del rp3_20a come from a WD+MS binary sys- 
tem evolved bv llvanova & Ta am (2004) with initial WD mass 
Mwd = 0.8 M . companion mass Md\ = 2.0 M G and orbital pe- 
riod P{ - Id (llvanova & Taaml [20041) . These parameters have 
been presented in Table 2 of PRWH08. For comparison, we set 
up our new Model ms_20c with the same initial binary param- 
eters (M WD = 0.8 M©, M d ,i = 2.0 M Q and P, = Id ) and carry 
out the fully detailed binary evolution calculation to construct 
the structure of companion star. Properties of Model ms_20c are 
listed in Table Q] The final orbital separation at the moment of 
the supernova explosion is 2.74 x 10 n cm, which is very close 
to the 2.68 x 10 1 1 cm of Model rp3_20a by PRWH08. Moreover, 
the final companion masses agree very well (1.22 M & for our 
Model ms_20c versus 1.17 M for Model rp3_20a). However, 
Fig. [8] shows that the ms_20c model has a larger radius com- 
pared to Model rp3_20a. Therefore, it is not surprising that the 
impact simulation for this model leads to a stripped mass of 
0.171 M (see TableQ]), while only 0.032 M were stripped in the 
rp3_20a model. Since our Model ms_20c has a larger radius, it 
has a less bound envelope that can be stripped away more easily. 



T— i— i— i— i— j- 
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Fig. 10. Fraction of the material stripped from the companion in 
velocity space for ms_22a model 5000 s after the supernova ex- 
plosion relative to the total contaminated supernova ejecta. The 
dashed line corresponds to the PRWH08 model. 



Moreover, the larger radius causes an extended interaction area 
that also leads to a slightly larger kick velocity of 64.32 km s 
compared with the 46.6 km s of Model rp3 J20a. Finally, in our 
model, a larger conical hole is created in the supernova debris be- 
hind the companion star. Without the effect of orbital separation, 
the more compact structure of the companion star significantly 
reduces the stripped mass. Thus, the degree of compactness of 
a companion star, especially the compactness of its outer layer 
is very important to determine the influence of the impact of the 
SN la ejecta on its binary companion star. 

What causes this difference of companion star structures 
between the models ms_20c and rp3_20a? PRWH08 computed 
their companion star models by constantly removing mass while 
they evolved a single main sequence star. They did not consider 
the detailed mass-transfer processes in a binary system. In their 
one-dimensional stellar evolution, the mass loss proceeds rather 
rapidly: the duration of the mass transfer is about a factor of 10 
less than the Kelvin-Helmholtz timescale of their stars. Although 
they mi micked the total mass lo ss and the total time of mass 
loss from llvanova & Taaml (120041) when constructing their com- 
panion star models, this does not lead to realistic companion 
star structures such as obtained from consistent binary evolu- 
tion. The result are more compact objects and the binary sys- 
tems of PRWH08 are characterized by very large values of the 
parameter, af/R2, since they used less evolved companion stars 
to start the mass accretion phase from, which were too compact 
to actually fill their Roche-lobe (see Fig. [6]). In contrast, RLOF is 
taken into account for the mass transfer of a binary system in our 
detailed binary evolution calculations, and we also consider the 
process o f WD accretio n by inc luding the optically thick wind 
model of lHachisu et al] dl999bl) . The structures of companion 
stars are properly adjusted due to the detailed modeling of mass 
transfer in the binary system. Consequently, our work presents 
an update to the PRWH08 study with more realistic companion 
star models. 
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Fig. 11. Relative fraction of material that belonged originally to the companion with respect to the total amount of material. The red 
and blue color correspond to companion and supernova material, respectively. 



In our SPH impact simulations, the stripped mass and the 
kick velocity depend sensitively on the ratio of the orbital sep- 
aration to the radius of companion star, More compact 
companion stars with a larger orbital separation lead to a signifi- 
cant reduction of the amount of stripped hydrogen-rich material. 
However, af/Ro is not the only parameter to determine the out- 
come of the supernova impact, the companion structure is also 
important (see Fig.|6]l. 



4.2. Distribution of the stripped material in velocity space 

The hydrogen-rich material stripped from the companion star is 
mixed into the supernova ejecta. Figure|9]shows the velocity dis- 
tribution of the hydrogen-rich material stripped from the com- 
panion star (left panel) and the supernova ejecta (right panel). 
Most of the stripped material is concentrated at velocities be- 
low ~ 800kms _1 , which is much slower than the typical veloc- 
ities of the superno va ejecta of ~ 10 000kms~ 1 dChugailll986t 
iMarietta et alJl200"oT) . placing it at the very center of the ejecta. 

In Fig.[l0]we plot the mass fraction of material stripped from 
the companion star in velocity space (all bound material has been 
cut out). We also show the result of PRWH08 (dashed line) for 
comparison. In the low velocity region, the mass fraction of the 
stripped material is very high and it dominates over the original 
supernova material. But this fraction sharply decreases as the 
velocity becomes larger than 10 3 kms~'. Clearly, the supernova 
ejecta dominate at high velocities. However, some stripped ma- 
terial is present also in the outer ejecta as already noted by MOO. 



They argued that the presence of stripped material at high ve- 
locities implies that traces of hydrogen from the companion are 
swept up in the oxygen and silicon layer of the supernova ejecta. 
They also argued that given the upper limits on the hydrogen 
abundance from SN la observations near maximum light, this 
may provide a criterion for discriminating between SN la pro- 
genitor scenarios (MOO). Figure[TT]illustrates the relative amount 
of material that originally belonged to the companion with re- 
spect to the total amount of material for different companion star 
models. This figure shows how the supernova ejecta are mixed 
with the material stripped from the surface of the companion 
star. Most of the stripped material, however, will only become 
visible at very late times when the ejecta are mostly transparent 
and it is possible to see very deep into the center of the ejecta. 

FigureQ~2]shows how the spatial distribution of stripped ma- 
terial evolves with time for model ms_20c. The left column 
shows both, the SN la and companion star at 1000 s, 3000 s and 
5000 s, respectively. The right column shows the material origi- 
nally belonging to the supernova only - all companion material 
has been cut out. We see that the supernova ejecta are signifi- 
cantly affected by the companion star when the supernova im- 
pacts the binary companion. The interaction creates a conical 
hole in the supernova debris with an opening angle of about 50° 
(see the right snapshot at 1000 s of Fig. \Y2\ . Comparing with 
PRWH08, our simulation shows a larger cone-like hole behind 
the companion star. This is not surprising, because at the same 
initial separation our companion star has a larger radius than the 
models used in PRWH08. Additionally, Fig. [12] shows that the 
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Fig. 12. Density distributions of all material (left panel) and the supernova ejecta (right panel) at 1000 s, 3000 s and 5000 s after the 
explosion for model ms_20c. The logarithm of density is color-coded. 



material is shocked at the companion star which leads to the for- 
mation of a bow shock. Ejecta passing through the shock are 
heated and com pressed into a thin shell, and their velocity vec- 
tor is redirected dKasenll2010h . 



4.3. Detection of hydrogen 

In our simulations, we find that the supernova impact strips off 
0.11 M to 0. 1 8 M Q of hydrogen-rich material from the compan- 
ion star (see Table 1). This is far more than t he most stringent 
upper limit of 0.01 M which iLeonar dl (l2007h derived from the 
non-detection of H ff emission in late time spectra. Therefore our 
results might challenge the SD scenario if the systems studied 
here are representative and the objects, from which the observa- 
tional constraints were derived, originate from SD progenitors. 

However, it is important to note that the model of 
iMattila et all d2005l) which was employed by Leonard! d2007l) 
to obtain the upper limit mentioned above is quite simple. 



Moreover, most of the stripped hydrogen in our models ends up 
at velocities below 10 3 kms~' so that it is confined to the inner- 
most part of the explosion ejecta which are usually rich in iron- 
group elements. Whether or not H a emission will be detectable 
under these conditions, is a highly complex question which can 
only be answered by performing sophisticated radiative transfer 
simulations on the abundance structure of our explosion models. 
This is beyond the scope of the current paper. 

Note also that there are some possibilities to reduce the 
amount of stripped hydrogen. In our det ailed binary evolution 
calculations, we adopt the prescription of lHachisu et al.l d!999al) 
for the mass growth of a CO WD by accretion of hydrogen-rich 
material from its companion. According to this model unpro- 
cessed matter is assumed to be lost from the binary system due to 
an optically thick wind if the mass accretion rate exceeds a criti- 
cal value of M C1 - WD = 5.3 x 10 7 (1.7/X- 1)(M WD -O.^Moyr" 1 
(here, X is the hyd rogen mass fraction and My p the mass of 
the accreting WD). lHachisu et"ai1 dl999al 120081) proposed that 



11 



Z. W. Liu et al.: 3D simulations of interaction between SN la ejecta and their MS companions 




Fig. 13. Density evolution of the remnant star in ms_22a model simulations. Only the bound material that originally belonged to the 
companion star are included. The logarithm of density is color-coded. 



this optically thick wind will strip off parts of the outer layer of 
the companion star. This "mass-stripping" effect is currently ne- 
glected in our simulations. However, it could reduce the amount 
of hydrogen in the SN ejecta since the mass-stripping reduces 
the companion size and therefore also the mass lost in the super- 
nova impact. 



Another possibility to reduce the amount of stripped hydro- 
gen arises if the MS comp anion in the binary has a helium-rich 
envelope (iPan et al.l (1201 Oh find the stripped helium mass to be 
consistent with observational constraints in this case). Such a 
system is possible if a massive primary undergoes a first RLOF 
during its red giant phase and the resulting helium star expe- 
riences a second RLO F episode during core helium burning 
dHachisuet al.lll999albh . 



Finally, we note tha t also the recently proposed "spin-up 
and spin-down" mo del dJusmaml 120111 15 i Stefano et alj|201 U 
lHachisu et al.ll2012h is likely to reduce hydrogen stripping. 



4.4. Surviving companion stars 

In the SD scenario of SNe la, the companion star survives the 
supernova explosion. Because the interaction with the SN la 
ejecta and its orbital velocity at the time of the explosion, its 
spatial velocity might distinguish it from stars in its neighbor- 
hood. In our simulation, the kick velocity reaches values from 
51 kms~' to 105 kms~', which is comparable to the orbital ve- 
locity of 96-281 km s _1 . Thus, the spatial velocity of the rem- 
nant star, v S patiai = -^ v kick + v orb' ran g es from 108 km s _1 (model 

ms_20 a) to 287kms _1 (model ms_30a). iRuiz-Lapuente et all 
d2004l) showed that Tycho G star has a spatial velocity of 
136 kms" 1 , which is located in the range of spatial velocities 
in our models. 

Furthermore, a surviving companion star will be strongly af- 
fected by the impact of an SN la and show distinguishing prop- 
erties. The identification of a surviving companion stars in his- 
torical supernova remnants is a promising method to test pro- 
genitor models of SNe la. Figure [13] shows the evolution of a 
remnant companion star up to 5000 s after the supernova explo- 
sion for model ms_22a. The remnant companion star is puffed 
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up due to the supernova impact and heating, and it has a large 
radius of ~ 4R G (see the last snapshot of Fig. [T3l ). Although 
the surviving companion star shrinks and relaxes to be almost 
spherical again 5000 s after the explosion, it is out of thermal 
equilibri um and its dens i ty and temperature are asymmetric (see 
Fig. [L3l) . lPodsiadlowskil d2003l) modeled the post-impact evolu- 
tion of the surviving companion star and showed that the star 
is able to completely recover thermal equilibrium ~ 500 yr af- 
ter the supernova explosion. He also pointed out that a surviving 
companion star may be significantly overluminous or underlumi- 
nous 10 3 - 1 4 yr after the explos ion relative to its pre-supernova 
luminosity (Podsiadlowski 2003). However, he did not simulate 
the dynamical interaction of the companion star with the super- 
nova ejecta. In our work, we followed the detailed interaction 
of SNe la and their companion stars. It is necessary to carry out 
a separate calculation to follow the post-impact evolution of a 
remnant star during its re-equilibration phase, which will con- 
strain the properties of the surviving companion star for searches 
in historic SN la remnants. This will be addressed in a forthcom- 
ing study. 

5. Summary and Conclusions 

Based on a range of more realistic companion star models than 
those used in previous work, we have investigated the impact 
of Type la supernova ejecta on their companion stars in WD 
+ MS binary systems using the SPH code GADGET3. As an 
initial model for our hydrodynamical impact simulations, all 
companion stars are constructed with Eggleton's stellar evolu- 
tion code incorporating the possibility of an instability of the 
accretion disk around the WD and including the prescription of 
lHachisu etai]dl999bl) for mass accretion onto the WD. We sum- 
marize the basic results and conclusions of our impact simula- 
tion as follows. 



1. 



For our binary systems, we always find a stripped mass larger 
than 0.1 Mq. This is in disagreement with the most recent ob- 
servational constraints on the detectio n of hydrogen in neb- 
ular spectra of SNe la dLeonar dl2007l) . 
Such large stripped masses cause a serious problem for the 
SD scenario of SNe la. However, prior mass- stripping from 
the companio n star by an optically thick wind dHachisu et al.l 
[l999a L 2008 ]) or a spin-up spin-down pha s e as proposed 



bv iJusthaml d201ll) : iDi Stefano etail d201ll) : lHachisu et all 
d2012l) might be able to explain the absence of hydrogen 
emission in nebular spectra of SNe la. 

3. For a given companion model, the dependence of stripped 
mass and kick velocity on the ratio of separation to the ra- 
dius of the companion, a^jRi, can be fitted by a power law. 
However, we do not find the same fit parameters as PRHW08 
when we put their models and ours together. This indicates 
that details of the structure of companion stars are important 
for the results of the supernova impact. 

4. The differences to previous works are attributed to a more 
realistic treatment of the binary evolution of our progenitor 
models. 

Further improvements to the binary evolution and observa- 
tional constraints for more SNe la are needed, and more detailed 
modeling (in particular of the radiative transfer in the nebular 
phase) is required to determine under which circumstances the 
hydrogen will be detectable in observations. 
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